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ABSTRACT 

We report the discovery of X-ray quasi-periodic oscillations (QPOs) at frequencies of 3-4 mHz from 
a transient ultraluminous X-ray source (ULX) X42.3+59 in M82. The QPOs are strong and broad 
and appear with weak or absent red noise, and are detected only in Chandra observations when the 
source is brighter than 10 40 ergs s _1 . The QPO behavior is similar to the type A-I QPOs found in 
XTE J1550— 564, which is a subclass of low frequency QPOs with properties in between type A and 
B. Therefore, we identify the QPOs in X42.3+59 as of type A or B, and rule out the possibility of 
type C. With this identification, the mass of the black hole in X42.3+59 can be inferred as in the 
range of 12,000-43,000 solar masses by scaling the QPO frequency to that of the type A/B QPOs 
in stellar mass black holes. Cool disk emission is detected in one Chandra observation, and the disk 
inner radius suggests a similar black hole mass range. Black holes of such a high mass are able to 
produce an energy output in a manner similar to X42.3+59 by accreting from the interstellar medium 
directly. 

Subject headings: black hole physics — accretion, accretion disks — X-rays: binaries — X-rays: 
individual (CXOM82 J095551.1+694045=X42.3+59) 



1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are variable X- 
ray sources which are not coincident with the nucleus of 
the galaxy and exhibit luminosities over the Eddington 
limit of a 20 M & black hole (3 x 10 39 ergs s" 1 ) assuming 
isotropic emission. They are too bright to be powered by 
accretion onto stellar mass black holes like Galactic black 
hole binaries (BHBs) and are thus candidates of interme- 
diate mass black holes (IMBHs). However, if the radia- 
tion is super-Eddington or beamed along the l i ne of sight, 
IMBH s are not required (Kin g et al.l l200ll iBegelmanl 
2002) . Modeling the X-ray spectra of ULXs could shed 
light on their masses. For example, the temperature and 
size of the accretion disk can be used to we igh the cen- 
tral compact object (TMakishim a et al1l2000[ ). However, 
identification and quantification of the disk emission in 
the energy spectrum, especially when it is the domi- 
nant component, has confront ed difficulties and been 
proved unreliable in a f ew ca ses (jGoncalves fc Soriall2006l : 
iFeng fc Kaaretl[2~007aL l2009h . In contrast, characteristic 
time scales of the X-ray emission, such as quasi-periodic 
oscillations (QPOs) and the frequency break in the power 
spectrum, could be used to determine the mass of com- 
pact objects via a model independent calibration (cf. 
iMcHardv et~aI1l2006t iCasella et al.ll2008t ). 

Low frequency QPOs from mHz to several tens Hz 
have been found in Galactic BHBs. On the basis of their 
properties like the coherence, amplitude, phase lag, and 
harmonic component, low frequency Q POs can be clas- 
sified into three types of A , B and C (fWiinand s et al.l 
[19991: IRemillard et aL|[200l . Type A/B QPOs appear 
in a relatively narrow frequency range at a few Hz, while 
type C QPOs vary in a wide frequency range in re- 
sponse to the spectral parameters of the source. For 
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example, the type C QPOs in GRS 191 5+105 are de- 
tecte d at frequencies of as low as 1 mHz (|Morgan et al.l 
Il997f h QPOs have bee n detected in two ULXs, and are 
argued to be of type C (jStrohmaver fc Mushotzkyll2"003l : 
iStrohmaver et al.l 120071 ) because they appear to vary in 
a manner similar to type C QPOs in Galactic BHBs. 
Scaling the compact object mass with the oscillation fre- 
quency indica^es_ttie_presmce of IMBHs in these two 
ULX s (|Feng fc Kaaretll2007ri IStrohmaver fc Mushotzkvl 
2009]). 

M82 is a starburst g alaxy at a distance of 3.63 Mpc 
(|Freedman et al.l I1994D . The X-ray source X42.3+59 
was identified as a transient ULX i n M82 from multi- 
ple X-ray observations with Chan dra dKaaret et al.ll2006l: 
IFeng fc Kaaretj2007bl : lkong et al.ll2007l ). IFeng fc Kaaret! 
(|2007bf ) argued that X42.3+59 is more likely to be an 
intermediate mass black hole than a stellar mass ob- 
ject accreting from a massive star according to its tran- 
sient nature (Kalo gera et al.l 120041) . It lies on the sky 
plane at 5" to another ULX X41.4+60, which at most 
times is the brightest source in M82, and thus Chandra 
is the o nly telescope able t o spatially resolve them in 
X-rays. iKaaret et al.l (|2006| ) reported significant timing 
noise near 1 mHz from X42.3+59 from one Chandra ob- 
servation. Here, using new Chandra data, we report the 
discovery of mHz QPOs in X42.3+59 and discuss their 
possible nature. 

2. OBSERVATIONS AND ANALYSIS 

We checked all archival Chandra and XMM-Newton 
observations of M82 to date, as well as three joint 
Chandra/XMM observations that we performed recently 
(PI: H. Feng), to search for timing noise of X42.3+59 
at low frequencies. Among these observations, QPO 
features at a few mHz are found in five observations 
(Chandra ObsIDs 6097/8190/10027 and XMM ObsIDs 
0112290201/0560590101), two out of which, 10027 and 
0560590101, were made simultaneously with an over lap 
of 17 ks. 
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TABLE 1 

Spectral and timing properties of X42.3+59. 



Instrument 


XMM 

(a) 


Chandra 
(b) 


Chandra 
(c) 


Chandra 
(d) 


XMM 
(e) 


Observation ID 
Observation date 
Exposure (ks) 
i?x (counts s — 1 ) 
Rs/Rt 


0112290201 
2001-05-06 
27 

2.59 

0.7 


6097 
2005-02-04 
53 
0.20 
1.0 


8190 
2007-06-02 
53 
0.31 
1.0 


10027 
2008-10-04 
19 
0.29 
1.0 


0560590101 
2008-10-03 
30 

3.27 

0.2 


Spectral Properties 


N H (10 22 cm- 2 ) 

r 

fx (10 -12 ergs -1 cm -2 s _1 ) 
L x (10 40 ergs s" 1 ) 
X 2 /dof 


2.35 - 2.47 


3.31 ±0.19 
1.44 ±0.09 
4.49 ±0.10 
1.13 ±0.04 
80.7/74 


3.22 ± 0.15 
1.31 ±0.07 
7.25 ± 0.12 
1.76 ±0.05 
87.1/74 


3.19 ±0.28 
1.33 ±0.13 
6.19 ±0.19 
1.51 ±0.08 
83.6/74 




QPO Properties 


Significance (<r) 
Frequency (mHz) 
FWHM (mHz) 
rms/mean (%) 


3.00 
3.98 ±0.18 
1.0 ±0.6 
5.0 ±0.9 


3.39 
2.89 ±0.35 
1.9 ±0.6 
10.2 ±1.6 


6.51 
3.58 ±0.21 
2.6 ± 0.5 
14.1 ± 1.1 


4.49 
2.77 ±0.33 
2.0 ±0.6 
17.4 ±2.0 


2.91 
2.97 ±0.33 
1.5 ±0.4 
14.4 ±2.6 



Note. — is the total count rate in the source extraction region, i?g/i?x is the fraction of source photons needed for the calculation 
of the QPO amplitude, TVjr is the absorption column density, T is the power-law photon index, /x is the observed flux in 1-8 keV, and Lx 
is the luminosity corrected for absorption in 1-8 keV. Errors are quoted at a confidence of 90% for spectral parameters and 68% (1<t) for 
timing parameters. 



TABLE 2 

Spectral fitting with a disk plus power-law model subject 
to absorption for observation 8190. 



Parameter 



Value 



T in 

Rln 

r 

^PL 

fx 
Lx 

disk fraction 
X 2 /dof 



(4.2 ± 0.4) x 10 22 cm' 

0.17 ±0.03 keV 

3.5±H x 10 4 km 

1.53 ±0.12 

(2.2 ±0.4) X 10~ 3 

(7.14 ±0.13) x 10- 12 

(2.7 ±0.7) x 10 40 

27% 

65.6/72 



Note. — Nu and T are the same as in Table [T] T- ln is the disk 
inner temperature, Ri n is the disk inner radius calculated assum- 
ing a face-on disk at a distance of 3.63 Mpc, TVpl is the power-law 
normalization in photons kcV -1 cm -2 s -1 at 1 keV, fx is the ab- 



sorbed flux in ergs 1 



Lx is the unabsorbed luminosity 



in ergs s 1 . The flux, luminosity, and the disk fraction are quoted 
in 1-8 keV. All errors are quoted at 90% confidence. 



For the three Chandra observations, the location of 
X42.3+59 on the charge coupled device (CCD) was more 
than 3.5' off the optical axis - where a point source 
spreads on multiple pixels and photon pileup for sources 
at the flux level of X42.3+59 is unimportant. Event 
lists were created using CIAO 4.1.2 with CALDB 4.1.3. 
Source energy spectra in the 1-8 keV band were extracted 
from the source region and grouped by a factor of 4 in 1-4 
keV, 8 in 4-6 keV, and 16 in 6-8 keV based on the origi- 
nal channels, with background subtracted from a nearby 
source-free region at the same location on the sky for the 
three observations. A power-law model subject to inter- 
stellar absorption provides adequate fits to the Chandra 
spectrum. The absorption column density is almost con- 
stant at about 3.3 x 10 22 cm' 2 . The source spectrum 
is hard, with a power- law photon index of about 1.3- 
1.4. The luminosity corrected for absorption varied by a 



factor 1.5 among the three observations. All spectral pa- 
rameters with 90% errors are listed in Table [T] We tried 
to add an additional multicolor disk component to the 
power-law spectrum. Pronounced improvement on the 
fits is only obtained for observation 8190; the disk com- 
ponent is favored at a significance of 4.1 a. The best- fit 
model parameters are listed in Table [5] with 90% errors 
quoted. Power spectra were calculated from lightcurves 
created using events from the same region for spectral 
analysis with a time step equal to the CCD frame time 
including the readout time, which is 0.44104 s for the 
three observations. Each lightcurve was divided into M 
segments with M = 25, 18, and 8 for observations 6097, 
8190, and 10027, respectively. A Fourier transform was 
performed for each segment and the individual power 
spectra were averaged and then geometrically binned by 
a factor of 1.2, 1.1, and 1.3, respectively for the three ob- 
servations. The binning factor is chosen for a reasonable 
frequency resolution around the QPOs and a sufficiently 
large binning factor needed for the fitting. We checked 
different energy bands and found that the most signifi- 
cant timing noise appears in the 1-8 keV range. 

For the two XMM observations, the lightcurves were 
produced using merged PN and MOS events in their 
common good time intervals (G TIs) from a half circu- 
lar region (region B in Fig. 2 of iFeng fc Kaaretl I2007U 
where the contamination from the other ULX X41.4+60 
is minimized) with a time step 10 times the PN frame 
time of about 0.73 s. Fourier transforms were performed 
in each GTI with 4096 points, and the power spectrum 
was obtained by averaging all individual ones. The num- 
ber of segments, M, is 8 and 9 and the geometric bin- 
ning factor is 1.1 and 1.2, respectively for 0112290201 
and 0560590101. For observation 0560590101, MOS1 
data were not included because they contain too many 
timing gaps. Energy spectra of X42.3+59 are not avail- 
able from XMM data due to considerable photon con- 
fusion with nearby sources. Spectral information for 
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0.001 0.010 0.100 

Frequency (Hz) 

Fig. 1. — X-ray Power spectra of X42.3+59 in the 1-8 keV energy 
range. The letters indicates the observation as shown in Table [T] 
The dotted line indicates the Gaussian plus cons tant model that fit 
the spectrum. The powers are normalized after [Leahy ct al., 111983 ) 
and the baselines of b, c, and d are shifted vertically for clarity by 
5, 10, and 15, respectively. 

XMM 0560590101 can be obtained from the simultane- 
ous Chandra observation 10027. For XMM 0112290201, 
which does not have a joint Chandra observation, the 
count rate of X4 2.3+59 has been resolv ed from surface 
brightness fitting (jFeng fc Kaaretll2007b| ) . There are four 
reliable spectral measurements of X42.3+59 with Chan- 
dra (the th ree ones reported here and observation 5644 
reported in lFeng fc KaaretJl2007bO . Taking into account 
the spectral variations of X42.3+59 with an absorption 
column density in the range of (3.2 — 3.4) x 10 22 cm~ 2 
and a power-law index in 1.3-1.5, the 1-8 keV luminosity 
of X42.3+59 in XMM 0112290201 can be estimated using 
PIMMS as in the range of (2.35 - 2.47) x 10 40 ergs s" 1 . 

Prominent timing noise above the Poisson level can be 
seen around a few mHz in the power spectra from these 
observations (Figure[T|), which we interpret as QPOs. We 
fit the power spectra in the whole frequency range with a 
a Gaussian plus constant model, which is overplotted on 
each spectrum with a dotted line. A Lorentzian also fits 
these peaks, but not as well as the Gaussian does. No 
obvious broadband red noise such as a (broken) power- 
law continuum is found near the QPO feature. We rule 
out the possibility that they are caused by instrumental 
oscillations like the dithering as follows. For the Chandra 
observations, photons from nearby regions do not show 
similar oscillations and simulations with MARX using 
the same aspect solution and bad pixel files confirms that 
no instrumentation signal is present at these frequencies. 



For the simultaneous Chandra and XMM observations, 
the same QPO feature is detected by both. The XMM 
power spectra are substantially contaminated by pho- 
tons from the nearby variable ULX X41.4+60, and the 
calculated QPO parameters may be affected. The pho- 
ton fraction of X42.3+59 in the extraction region is esti- 
mated and taken into account when calculating the rms 
amplitude of QPOs from the XMM data. All QPO pa- 
rameters are listed in Table [I] and the QPO significance 
is calculated based on the x 2 variation on the normal- 
ization of the QPO component. We checked the power 
spectra at even lower frequencies for the three Chandra 
observations, and found possible red noise only in 6097 
appearing at frequencies well below the QPO frequency. 
For observation 8190, the power spectrum between 0.01- 
1 mHz is dominated by the white noise. For observation 
10027, which is shorter than the other two, the powers 
between 0.1-1 mHz are also white noise dominated. 

There are several other Chandra Advanced CCD 
Imaging Spectrometer ( ACIS) observations i n which 
X42.3+59 was active (see lFeng fc Kaareti[2007bl ) but the 
mHz QPOs were not detected. Observation 5644 has 
an exposure of 70 ks and the source luminosity was 
9 x 10 39 ergs s _1 in 1-8 keV. The number of photons 
from X42.3+59 detected in 5644 is as many as in 6097, 
and twice that in 10027; the nondetection of QPOs is not 
due to low statistics. Observation 2933 and 6361 both 
have an exposure of 18 ks and the source luminosity was 
4 x 10 39 and 8 x 10 39 ergs s _1 , respectively. Another 
four ACIS observations (378, 379, 380-1 and 380-2) are 
shorter than 10 ks, and the source was slightly dimmer 
than in observation 2933. The mHz QPOs were not de- 
tected in any observation when the source was dim, below 
- 10 40 ergs s _1 . 

3. DISCUSSION 

Low frequency QPOs in stellar mass black hole X-ray 
binaries have a diverse population distributed in a wide 
frequency range. Therefore, inference of the compact ob- 
ject mass via the QPO frequency should be done between 
QPOs of the same type. The QPO coherence parameter 
Q, defined as the ratio of the central frequency to the 
full width at half maximum (FWHM), is around 1.4-1.5 
derived from the three Chandra observations and nearly 
4 from the XMM observation 0112290201. The rms am- 
plitude of the QPOs in the three Chandra observations 
is high, from 10% to 17%, but is relatively low in the 
XMM observation 0112290201 of 5%. Wc note that the 
Chandra results should be more reliable due to less con- 
tamination. As mentioned above, possible red noise is 
only found in Chandra 6097, and the power spectrum is 
totally dominated by the white noise in the other two ob- 
servations at frequencies down to 0.1 or even 0.01 mHz. 
Therefore, the QPOs in X42.3+59 can be summarized 
as being broad and strong, and occurring with absent or 
weak red noise continuum. 

The QPOs in X42.3+59 are similar to four QPOs de- 
tected in XTE J1550-564 on 1999 Mar 18 and 21 and 
Apr 02 and 03, respectively, which also have a low Q 
of about 1.0, an rms amplitude of about 10%, and ac - 
companied with a weak red noise (jRemillard et all20 02) . 
Another similarity is that three out of the four QPOs in 
XTE J1550— 564 do not exhibit a Lorentz profile. These 
four QPOs were classified as 'possible' type A QPOs by 
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IRemillard et al.l (12002H. o r type A-I (a subclass of type 
A) bv lHoman et al.l (12001|) to distin guish them from weak 
type A QPOs. iCasella et all (|2005l ) argued that type A- 
I QPOs should be classified as type B based on their 
frequency and amplitude. However, type B QPOs are 
usually narro w, with a Q value o f > 6 concluded from 
a few BHBs (|Casella et al.ll200l . Therefore, the type 
A-I could be a special subclass amid type A and type B. 
In spite of the controversy over the classification of type 
A-I, we identify the QPOs found in X42.3+59 as of type 
A-I based on their similar profile and the weak/absent 
red noise, or conservatively, of type A/B according to 
the discussion above. As noted above, the mHz QPOs 
were detected only when the source was luminous, above 
~10 40 ergs s -1 . Type A/B QPOs show identical be- 
havior that they appear only w hen the source flux is 
high enough IjCasella et al.l l2004). This further supports 
the identification of the mHz QPOs as of type A/B. By 
any means, the QPOs in X42.3+59 are not analogous to 
type C, which are narrow, strong, and accompanied by 
a strong red noise continuum; otherwise, the following 
discussion would be different as type C QPOs appear at 
a frequency range significantly wider than others. 

The four type A-I QPOs in XTE J1550-564 have a 
central frequency varying between 7.1 to 10.3 Hz; all 
type A QPOs in XTE J1550-564 are found in a fre- 
quency range of 4.88-10.3 Hz, which covers the frequency 
range of type B QPOs (4.94-6.1 Hz) in the same source 
(|Remillard et all 12002( 1 . The QPOs in X42.3+59 have 
a central frequency between 2.77-3.98 mHz. Comparing 
with the type A-I QPOs in XTE J1550-564, the upper 
and lower bounds of the frequency range have exactly the 
same ratio but scaled by a factor of about 2600. Based 
on the assumption of type A-I QPOs, the compact ob- 
ject mass of X42.3+59 can be inferred as 25,000-30,00 
M Q after a dopting a mass of 9.68-11.58 M for XTE 
J1550-564 ((O rosz more conservative es- 

timate taking into account all type A/B QPOs leads to 
a black hole mass of 12,000-43,000 M s in X42.3+59, in- 
ferred from the possibly smallest and largest ratios be- 
tween QPO frequencies from the two sources. 

Type A/B QPOs are detected in several other BHBs 
besides XTE J1550-564, including GX 339-4 XTE 
J1859+226, and GS 1124-684 f see ICasella eTall 120051 
and references therein). In these sources, the type A/B 
QPOs show up in a similar frequency range, e.g. 4.5-7.49 
Hz in GX 339-4 (IBelloni et al.l 20051) and 4.45-8.42 Hz in 
XTE J1859+226 (Cas ella et al.ll2004[) . Type B QPOs are 
also detected in GRS 1915+105 varying in the frequency 
range of 2.44-6.84 Hz (|Soleri et al.ll2008D . They are broad 
but shown above a strong red noise component, par- 
tially similar to the QPOs in X42.3+59. These sources 
also contains a similar massive black hole as in XTE 
J1550— 564, but their masses are not as well con strained 
as in the latter (jRemillard fc McClintockl 12006(1 . There- 
fore, inferring the mass of X42.3+59 via type A/B QPOs 
from these sources will lead to consistent but looser con- 
straints. We note that all above estimates of the mass 
depend on a hypothesis that the frequency of type A/B 
QPOs is linearly scaled with the black hole mass, which, 
however, has never been tested. More work regarding 
the mass dependence of type A/B QPOs is needed. 

In GX 339-4, all type A/B QPOs were detected 
when the source was in the soft intermediate state 



(|Belloni et al.l 120051 ) , also known as the steep pq wer- 
law state defined by IRemillard fc McClintockl (|2006H . In 
XTE J1550-564, the four type A- I QPOs were found in 
the steep power-law state as well (|Remillard et al.ll2002l : 
iSobczak et al.l 12000( 1 . However, the X-ray spectrum of 
X42.3+59 is hard, seemingly inconsistent with the steep 
power-law state. Disk emission is detected in observation 
8190, with a moderate fraction of 27% in 1-8 keV. This 
disk component will contribute 40% to the total emission 
below 1 MeV if the power-law component is not cut off. 
In the hard state, the disk contributes less than 25% in 
2-20 keV for stellar mass black holes. For more massive 
black holes, though there is no defined energy range, a 
disk fraction of 27%-40% seems too high for being clas- 
sified as in the hard state. A varying disk component 
and an energetic power-law component suggest that the 
source may be in a state similar to the steep power-law 
state in stellar mass black holes. Interestingly, the best- 
fit disk inner radius is consistent with the gravitational 
radius of a black hole of 11,000-44,000 M Q , which is in 
agreement with the mass estimated from the QPO fre- 
quency. 

X42 .3+59 lies about 5" fr om the kinematic center of 
M82 (|Weliachew et all 11984ft on the sky plane. Follow- 
ing the discussion in lKaaret et al.l (|2001l ) and references 
therein, an upper limit of the compact object mass can 
be placed of about 3 x 10 4 M & such that the black hole 
will not fall onto the nucleus of the galaxy due to dy- 
namical friction if the black hole was born at the same 
time of the galaxy. This upper limit is consistent with 
our mass estimate via QPOs. 

X42.3+59 is coincident with an extended radio source 
known as 42.21+59.2, which i s believed to be an H II 
region (|McDonald et al.ll200l: iKaaret et al] [20061) asso- 
ciated with a star cluster seen from infrared (jKong et al.l 
120071 ). Therefore, it could be ruled out that the source 
is a background AGN. If the disk component is true, its 
high bolometric luminosity of 1.2 x 10 41 ergs s _1 and low 
temperature make it a good candidate for the ionization 
source of the H II region. 

The outburst luminosity of the source is about 
10 40 ergs s . In the total 20 observatio ns of the source 
(16 reported in iFeng fc Kaaret) l2007bl plus Chandra 
8190, 10027, 10025, and 10026), X42.3+59 was found ac- 
tive in 12 of them, indicative of a duty cycle e = 0.6. As- 
suming an accretion efficiency rj = 0.1, the mean mass ac- 
cretion rate of the source is found to be M = eLx/i]c 2 
7 x 10 19 g s _1 . If the source is accreting from the inter- 
stellar medium directly, the mass accretion rate in g s" 1 
can be written as 

m = 1.4x10" ("X ( j \ ( c - y 

\MqJ Vl0- 24 gcm-V U0kms-V , 

where M is the black hole mass, p is density and c s i s 
the speed of sound of the medium dFrank et al.l 12002). 
Adopting canonical values of the density and sound speed 
shown in the equation, a black hole mass of 2 x 10 4 M Q 
is required to power the X-ray luminosity without a com- 
panion star. This means that if the mass estimated from 
the QPOs are reliable, the source does not have to be in 
a close binary system, but is massive enough to a ccrete 
from the interstellar medium directly. iPellegriml (|2005H 
suggests that the Bondi rate may overestimate the true 
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accretion rate for AGNs by 1-2 orders of magnitude. In 
this case, a corresponding higher black hole mass would 
be required. Conversely, if the ULX indeed lies within a 
star cluster, the gas density could be substantially higher, 
which would decrease the required black hole mass. Any- 
way, one should be cautious of the high uncertainty on 
such an estimate. 
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